Finite Temperature Magnetism in Fractional Quantum Hall Systems: Composite 

Fermion Hartree-Fock and Beyond 
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Using the Hamiltonian formulation of Composite Fermions developed recently, the temperature 
dependence of the spin polarization is computed for the translationally invariant fractional quan- 
tum Hall states aX v = 1/3 and u — 2/5 in two steps. In the first step, the effect of particle-hole 
excitations on the spin polarization is computed in a Composite Fermion Hartree-Fock approxima- 
tion. The computed magnetization for i/ = 1/3 lies above the experimental results for intermediate 
temperatures indicating the importance of long wavelength spin fluctuations which are not correctly 
treated in Hartree-Fock. In the second step, spin fluctuations beyond Hartree-Fock are included 
for v = 1/3 by mapping the problem on to the coarse-grained continuum quantum ferromagnet. 
The parameters of the effective continuum quantum ferromagnet description are extracted from the 
preceding Hartree-Fock analysis. After the inclusion of spin fluctuations in a large- Ai' approach, the 
results for the finite-temperature spin polarization are in quite good agreement with the experiments. 
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I. INTRODUCTION 

The fractional quantum Hall (FQH) efFectfl has in- 
troduced us to new, highly correlated, incompressible 
statesH of electrons in high magnetic fields. A unified un- 
derstanding of all fractions v = p/{2sp+ U was achieved 
by the Composite Fermion picture of Jainlj, in which the 
electrons are dressed by 2s units of statistical flux to 
form Composite Fermions (CFs). At a mean field level, 
the CFs see a reduced field B* = B/{2sp+ 1), in which 
they fill p CF-Landau levels (CF-LLs), and exhibit the 
integer quantum Hall effect. 

Due to the small g factor of electrons riip. GaAs, spins 
may not be fully polarized in FQH statesQEI. Transitions 
between singlet, partially polarized, and fully polarized 
states (based on gap measurements) have been observed 
for a number of fillingsQ^ pathich can be understood in 
terms of CF's with a spinl3'BE3. The transitions happen 
when an unoccupied CF-LL of one spin crosses the occu- 
pied CF-LL of the opposite spin. 

While these low temperature measurements are in sat- 
isfactory agreement with the ground states predicted in 
the CF picturglj, in order to understand the tempera- 
ture dependence of the polarization P{T) one has to con- 
sider all excited states as well. Detailed measurements of 
P{T) for pk9-i/ = 1/3 state have recently appeared in the 
literaturey'tHl. It is well-known that the v — 1/3 state 
is spontaneously polarized at T = 0, even when the Zee- 
man coupling Ez ^giiBtot is zero. In this it is analogous 
to the v = l|-statfiE3, which has been |ej;tensively studied 
theoreticallyliJ~E3 and experimentallyO. There are, how- 
ever, significant differences between the two cases at finite 
T. 

In a recent paper, MacDonald and PalaciosEfl identi- 
fied a key qualitative feature that makes v = 1/3 very 
different from ly — 1. In the 1^=1 case the particle- hole 
excitations are very high in energy compared to Ez , and 



are frozen out at all low temperatures of interest. Con- 
sequently, the T dependence of P comes mainly from 
spin wave excitations and their interactions. This is the 
reason why long- wavelength effective theories such as the 
continuum quantum ferromagneiliJ approach are success- 
ful. However, for v = 1/3, particle- hole excitations are on 
the same scale of energy as Ez , and cannot be ignored at 
any T. MacDonald and. Palacios use a simplified model to 
illustrate this featureEfl, but the model is not sufficiently 
detailed to enable a calculation of P{T) for a realistic 
sample, and is difficult to extend to non-Laughlin frac- 
tions. 

The goal of this paper is to describe a general an- 
alytical method for approximately computing P{T) for 
an arbitrary principal fraction for realistic samples. To 
this end, an approximate hamiltonian formalism in which 
Conippsite Fermion variables explicitly appear will be 
usedllj. This hamiltonian approach is based on the field 
theoretic idea of aitaphing flux to electrons by using a 
Chern-Simons fiel(£j"E3. The CF-Hamiltonian approach 
has many features suited to the computation of the phys- 
ical properties of fractional Hall systems. The CFs see 
the effective field, and fill CF-Landau levels (CF-LLs). In 
the principal fractions v — p/ (2p-|- 1) with which we will 
be concerned, the effective field is B* = B/{2p+ 1), and 
the CFs fiU p CF-LLs. The energies of the CF-LLs are 
controlled entirely by interactionstJ, which is a correct 
feature of the physics of electrons in the lowest Landau 
level. Note that this is not a theory in which Composite 
Fermions are free. On the contrary, the theory is fully in- 
teracting, with both the kinetic energies and the residual 
interactions of the CFs being controlled by the electron- 
electron interaction. Finally, the nonperturbative charge 
and dipoleJBoment of the excitations appear explicitly in 
the theoryt£l. The fact that all these nontrivial features 
are built into the theory raises the expectation that very 
simple approximate treatments of this Hamiltonian the- 
ory (such as Hartree-Fock) will suffice to produce reason- 
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able numbers for physical quantities. This exoertation 
is indeed borne out by explicit calculation£3LJ. The 
hamiltonian formalism will be presented briefly in Sec- 
tion II. 

The Hamiltonian approach is particularly suited to 
the computation of finite-temperature properties in frac- 
tional Hall systems. In the Composite Fermion Hartree- 
Fock (CFHF) approximation, one self-consistently finds 
the single-particle energies and occupations of the various 
CF-LLs at any finite temperature. Since the energies of 
the states are controlled entirely by interactions and oc- 
cupations, these energies will be temperature- dependent. 
This is a familiar feature in other interacting many-body 
systems, such as the BCS superconductor, where the 
single-particle gap is a collective effect, and depends on 
temperature. Thus the Hamiltonian theory gives us valu- 
able information on the evolution of the collective state as 
a function of temperature, which is then used in the map- 
ping to the effective theory. Once the occupations of the 
single particle CF-LLs have been determined, the polar- 
ization is computed simply as the difference of the total 
"f-spin and J,-spin occupations. This yields the CFHF pre- 
diction for the temperature dependence of the spin mag- 
netization. This procedure and the results for = 1/3 
and V = 2/h are described in Section III. j4J3rief descrip- 
tion of this work has appeared previouslyE3. 

It turns out that for the spontaneously polarized v = 
1/3 state the CFHF prediction is higher than the ex- 
perimental values for all temperatures, with the discrep- 
ancy being considerable for intermediate temperatures, 
as shown in Figure 1. To put it in the proper context, 
the agreement between CFHF and experiment is consid- 
erably better than for v — \, where HE_does a very poor 
job of predicting the spin polarizatiorilj. 
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FIG. 1. Polarization versus T for = 1/3. The citeJes are 
the data from the lOW sample of Khandelwal et aiiJ. The 
dashed line is the prediction from the CFHF theory for a 
thickness parameter of A = 1.5L The sohd and dot-dashed 
lines refer to the theoretical prediction including spin fluctu- 
ations. 

However, the discrepancy is nonetheless there, and is 
presumably the result of spin fluctuations which are not 



treated correctly in the CFHF approximation. One can 
imagine integrating out the fermions to obtain an effec- 
tive theory that has low energy spin degrees of freedom. 
One is then led to map the low energy physics, on to 
the continuum quantum ferromagnet (CQFM)EJ. The 
CQFM has two free parameters, the magnetization per 
unit volume Mq, and the spin stiffness ps. In the tradi- 
tional CQFM theory these are temperature-independent 
parameters. However, since the theory is fermionic at the 
microscopic level, with the fermionic energy levels and oc- 
cupations being temperature dependent, one should ex- 
pect Mq and ps to acquire a T dependence in the effective 
theory. It turns out that these parameters can be eas- 
ily extracted from the CFHF treatment of Section HI. 
Armed with this information, we proceed to include spin 
fluctuations ia. a large-N approach as described by Read 
and SachdevEj. The results are in quite good agreement 
with the experimental data over the whole range of tem- 
perature, as shown in Figure 1 for the lOW sample of 
Khandelwal et aO. Thei-s ame comparison is shown for 
the data of Melinte et aO in Figure 2. 
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FIG. 2. The theoretical prediction in the CFHF approx- 
imation (dashed line) and including spin fluctuations (solid 
and dot-dashed lines) compared torthe Knight shift data from 
the M242 sample of Melinte et ali3. A saturation value of 
19kHz for the Knight shift has been used, and the thickness 
paramter is assumed to be A = 1.5/. 

The details of the figures will be explained in Sections 
III and IV. The mapping to the effective theory and the 
subsequent calculations are described in Section IV. We 
end with some conclusions, caveats, and open questions. 



II. HAMILTONIAN FORMALISM 

Since detailed descriptions of |the Hamiltonian theory 
of CFs have appeared elsewheralj, we will restrict our- 
selves to a summary of the essential features of this for- 
malism. 

The ^Composite Fermion picture was originally used 
by JainEl to generate electronic wave functions with good 
correlations. These wave functions have excellent over- 
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lap with the exact wave functions for finite systems, and 
encode all the right physics. In order to compute dynam- 
ical response functions, it is desirable to have an operator 
or field theoretic formulation of CFs. The fundamental 
property that CFs carry statistical flux can be imple- 
mented by a Chern-Simons(CS) transformation, which 
performs flux attachment via the CS gauge field to-pb^ 
tain a field jtljeoretic description with either bosonsEHTo 
or fermionaSa. These theories have provided us with a 
link between the microscopic formulation of the problem 
and experpaeijit, both for incompressible and compress- 
ible statesEaEa. 

Recently R. Shankar and the present author developed 
a hanuljtonian CS theory for the fractional quantuniJIall 
stateslij. Inspired by the work of Bohm and PineS on 
the 3D electron gas, the Hilbert space was enlarged to in- 
troduce n high-energy magnetoplasmons degrees of free- 
dom, (n also being the number of electrons) at the same 
time imposing an equal number of constraints on physical 
states. However, the fermions still had the bare mass (re- 
call that in the lowest Landau level, the electrons should 
lose all memory of the bare mass, and acquire an effective 
mass controlled by interactions), and the frequency of the 
magnetoplasmons was incorrect. Hence a final canonical 
transformation was employed to decouple the fermions 
from the oscillators in the infrared limit. 

The final fermions are called the Composite Fermions 
for the following reasons. Firstly, the final fermions have 
no dispersion in the absence of interactions and acquire 
an effective mass dependent on interactions alone. Next, 
the final canonical transformation assigns to each fermion 
the magneiic-inoment e /2m as mandated by very general 
argument£3o. The central result of the formalism is 
the formula for the electronic charge density, which takes 
the following form, separable into high- and low-energy 
piecesEj, at small q: 



pM = -j=J-^{A{q) + A\-q)) 
VStt V 2p + 1 



2p+l 



(1) 



where A, A^ refer to the annihilation and creation oper- 
ators of the magnetoplasmon oscillators, I — l/VeB is 
the magnetic length, and Hj = Pj + eA*{rj) is the ve- 
locity operator p£the CFs. The oscillator piece saturates 
Kohn's theoremEj. The rest, to be called p, is obtained by 
adding to the canonically transformed electronic charge 
density a particular multiple of the constraintlij (in the 
physical subspace, one can add any multiple of the con- 
straint without physical consequences, but we wish to 
work in the full space). It has some very useful proper- 
ties in the full space: 

• p satisfies the magnetic translation algebra 
(MTA)Ej to lowest leading order. Since this is the 
algebra of the electron density in the lowest Landau 
level (LLL), the LLL projection has been correctly 
carried out in the infrared. 



• Note that p is a sum of a monopole with charge 
e* — e/ (2p-|- 1), which is the charge associated with 
the CF, and a dipole piece which alone surviveSpat 
v = 1/2 and has the value proposed by ReadEJ. 
(A number of recent-coDstructions have emphasized 
this dipolar aspecttZrc3) . 

• We also find that that asq^O all transition matrix 
elements of p from the HF ground state vanish at 
least as q^. 

The final property is an essential property of physical 
charge density matrix elemerita from incompressible liq- 
uid ground states in the LLLEj. It is easy to see that if 
one intends to use the Hartree-Fock approximation ignor- 
ing constraints, these properties of p are essential. They 
make it plausible that p does not suffer vertex correc- 
tions. 

The Hamiltonian of the low-energy sector (dropping 
the magnetic moment term) is 



H 



v{q)p{-q)p{q)- 



(2) 



where v(q) is the electron-electron interaction. Real 
samples have a finite thickness A of the same order as 
I, so that ttje Coulomb interaction is cutoff at, large 
wavevectorscll. We will model this interaction byE3 



v{q) 



2^e2 



■ exp —Xq 



(3) 



Finally, the constraint will be approximately imple- 
mented by cutting off the number of CF-LLs to maintain 
the correct number of electronic states. For v — pj (2p-(-l) 
this means keeping 2p -I- 1 CF-LLs. More details of the 
formalism can be found in refs. [19,29]. 



III. THE COMPOSITE FERMION 
HARTREE-FOCK APPROXIMATION 

The CFHF approximation has been applied to the 
above Hamiltonian, and reasonable success has been ob- 
tained in computirig,. various physical quantities in the 
gapped fractionsEZrHj, including a very recent calcula- 
tion of the temperature-dependent polarization P(r) for 
the compressible half-filled LLcj. 

Before one employs the CFHF approximation, one 
needs to express the Hamiltonian as an operator acting 
on a set of states. Since the CFs see the effective field 
S*, it is natural to represent the Hilbert space as Slater 
determinants of single-particle CF-Landau level states in 
the effective field. The wave functions of these states in 
the Landau gauge are 



(4) 



Here I* — l\J2p -I- 1 is the effective magnetic length, L is 
the linear dimension of the system, X ~ 2'k{1*Y j / L is the 
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degneracy index, and j = 0, 1, L^/(27r(r)^), and (/>„ 
is the n*'' normalized harmonic oscillator eigenfunction. 
The index n is the CF-LL index. 

The density is a one-body operator, and can be ex- 
pressed in this basis (after the spin labels have been in- 
cluded) as 



Pig) 



E 

cr,X{ni} 



<T,ni .X- 



Mllr. cr,n2,X + 



(5) 



where a =t, j is the spin index, and da-^n,x destroys a 
CF in the single-particle state labelled by cr, ri,X, and 
p„j„2((z) is a matrix element given by 



M 



(-1)" 



Pnin2 \HJ — 2p+l 



n > — n < 
1 




n("i-"2)(-£'i)">- 



(n< + l)L::>;r) (6) 



where (n^) is the lesser (greater) of ni,n2, Og is 
the angle of the vector q measured from the x-axis, 
y = {ql*Y /2 is the argument of the Laguerre polyno- 
mials L^. 

The Hamiltonian is now a four-fermi operator. It can 
easily be shown that the siagle-particle states defined 
above form a good HF basiscJ. In the CFHF approxima- 
tion, one reduces the four-fermi Hamiltonian to a two- 
fermi operator by taking averages according to the rules 



< GS\dld^,\GS >= 5^^,Nf{v) 
< GS\d,dl,\GS S,,,{1 - Nf{i^)) 

where ;/ is a shorthand for all the state labels, 
results in the HF Hamiltonian 



HF 



cy,n,X 



(7) 
(8) 

This 



(9) 



where the HF single-particle energy is 
e(CT, n) = 

± ^ + U ErJl - ^Np{a, m))|p_(g)p (10) 

in which the Zeeman energy has been added {Ez = 
g*fiBB). 

Finding the energies at T = is quite simple, since the 
occupations of the CF-LLs can only be or 1. For ex- 
ample, for the u = 1/3 state, Np{1,0) = 1 and all other 
occupations are zero, and for the v — Ijh singlet state 
Af_F(t,0) = A^F(i,0) = 1, and all other occupations are 
zero. For nonzero T one has to carry out a self-consistent 
procedure. First one chooses trial values for the energies 
(say the T = values) and a trial value for chemical 
potential /i (say halfway between the lowest unoccupied 
and highest occupied CF-LL). Then one assigns the oc- 
cupations of the single-particle levels according to Fermi 
occupation function 



1 



1 -I- exp (e(cr, n) — [i)lT 



(11) 



with the trial value of [i. Then one recomputes the HF 
energies using equation ([l0|). Note that the structure of 
degenerate CF-LLs remains intact, and only the energies 
and occupations change. Since the filling has to remain 
fixed, the chemical potential will change when the ener- 
gies change. One then recomputes the chemical potential 
as the root of the equation 



(12) 



for the principal fraction p/ (2p-|- 1) and iterates the whole 
process until self-consistency is acheived. Finally, the 
spin polarization is given by 



P = ^(7Vj.(t,n)-iVj.(i,n))/p 



(13) 



From the above procedure it is clear that only single- 
particle excitations have been taken into account in ob- 
taining P(T) so far, and spin fluctuations have been ex- 
plicitly ignored. If it so happens that for the system 
under consideration the effects of spin fluctuations are 
small, then this approximation should be accurate, oth- 
erwise not. 

Let us proceed to compare the CFHF results to ex- 
periments. pWe first consider the lOW sample of Khan- 
delwal et aO. The sample parameters are = 9. GIT, 
and Btot = 12T. This implies that the Coulomb energy 
scale is Ec = jel « IGOiiT and the Zeeman energy is 
Ez = 0.0175£'c. We wih use a value of A = \.hl for the 
thickness parameter for illustrative purposes through- 
out this-. paper. This value ought to be in the physical 
regimeE3 for most samples, and also approximately agrees 
with that extracted from an analysis of the compressible 
statescj. It should be emphasized that the CFHF anal- 
ysis, and the mapping to the effective spin theory that 
follows, can he performed for any potential v{q). 

Figure 1 (in the Introduction) shows the HF predic- 
tion from our theory for A = 1.5/ compared to the ex- 
perimental data. The agreement is good at very low 
[T < IK) and very high (T > 8K) temperatures. How- 
ever, at intermediate T, there is a big discrepancy be- 
tween the CFHF prediction and the data. This indicates 
that effects not treated correctly in HF, notably long 
wavelength thermal spin fluctuations, are important in 
this intermediate regime of T. Nonetheless, to put the 
result in context, one should note that the CFHF pre- 
diction agrees much better with the data than the cor- 
responding HF prediction for v ~ 1 (see, for example, 
ref[15] for a comparison of the different predictions for 
V = 1). This is because at = 1 particle- hole excitations 
are completely unimportant at all temperatures of inter- 
est, while thermal spin fluctuations dominate. Since long 
wavelength spin fluctuations are treated very poorly in 
HF the agreement is bad. However, at = 1/3, particle- 
hole excitations play a major role in reducing P{T) for 
T > 5K. In the next section we will see how to incor- 
porate spin fluctuations into our calculation, resulting in 
much better agreement with the data. 
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FIG. 3. Knight shift versus T for v = 1/3. The circles are 
the data from the M242 sample of Melinte et ao. The lines 
are the predictions from the CFHF theory for A = 1.5Z, assum- 
ing two different values for the Knight shift that corresponds 
to saturated polarization. 



the following interesting question: Why is HF so good 
in this case relative to the case of = 1? To answer 
this question let us turn to the spin wave dispersions. 
The spin wave is a collective spin-flip excitation, and at 
wave vector q corresponds a plane wave state in which a 
majority spin quasihole and a minority spin quasiparticle 
are at a separation of q{r*)'^. 
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Figure 3 shows the same type xtf comparison for the 
Knight shift data of Melinte et aOfor their M242 sam- 
ple. Here the sample parameters are Btot — B± = 17T. 
This implies that Ec « 210K, Ez = O.OlSGSc ~ AK. 
Once again a value of A = 1.5Z has been used in the CFHF 
calculation. In order to translate the Knight shifts into 
polarization numbers or vice versa, one needs to deter- 
mine the saturation value of the Knight shift, which pre- 
sumably corresponds to a polarization of P = 1. There 
is a lot of scatter in the data at low T, due to the very 
long times needed to measure the_Knight shift, and the 
error bars are also large at low Tli3. There is thus some 
uncertainty in determining the Knight shift correspond- 
ing to P = 1. In any reasonable theory one expects to 
find that P = 1 for T ^ Ez, and expects to see this 
saturated value of P up to about T — 0.5Ez or so. 

Based on these considerations two values of the satu- 
ration Knight shift Ks^p=i, 21kHz and 19kHz have been 
used to fit the data here, both of which lie within the er- 
ror bars of the low T dataHl. One possibility that can ex- 
plain this spread is that spin-reversed quasiparticles are 
present in the ground state due to disorder, which qaii 
bring down the "saturated" value of the Knight shifted. 
The 21kHz value was used by Melinte et al in a phe- 
nomenological tanh(A/4fcBr) fit to obtain A = l.lEz- 
However, one must note that the fit for Ks^p=i = 19kHz 
seems slightly better, since then the experimental satu- 
ration region is about O.bEz- The agreement between 
theory and experiment for this value of K^^p^i are also 
better than for Ks^p=i — 21kHz. The 19kHz value will 
be used in the mapping to the effective theory in the 
next section. Overall the agreepajent is slightly worse than 
for the Khandelwal et al datallil, but leads to the same 
conclusion: It is quite important to treat thermal spin 
fluctuations correctly at intermediate temperatures for 
ly ^ 1/3. 

Before we proceed to approximately incorporate spin 
fluctuations into the theoretical prediction, let us address 



FIG. 4. Spin wave energy dispersions in units of the 
Coulomb energy Ec for Ez — 0.0175 and A — 1.51, with 
q plotted in units of . As can be seen, the scale of the 
spin-reversed gap is the same as Ez. 



The g — > limit is required to be Ez by Larmor's 
theorem, while in the q —^ oo limit the particle and 
hole become infinitely separated, so that the energy of 
the excitation is the spin-reversed particle-hole gap Asp. 
The dispersion of these excitations can be computed for 
J/ = 1/3 in the manner described in ref . [29] , and is shown 
in Figure 4 for A = 1.5^ for Ez = O.OllMJc and P = 0. 
Figure 4 explicitly illustrates the featuretS that the spin- 
flip particle-hole excitations are at the same energy scale 
as Ez- 

How does this evolve as temperature increases? Recall 
that all the energy splittings in the CF-Hamiltonian for- 
malism come from interactions, and as the occupations 
of the states change with T so do their HF energies. As 
T increases the occupations of the minority spin levels 
increase while that of the lowest majority spin level de- 
creases. This means that the exchange splitting between 
the minority and majority spin levels decreases, as can 
be seen from equation (|o|) . It is clear that as T becomes 
very large the occupations of all the levels should tend 
to become the same, and therefore Asp^ should tend to- 
wards Ez- Since the spin wave dispersion has to be Ez 
for very small q, and Asp for very large q, this implies 
that the spin wave dispersion becomes increasingly flat 
at T increases. We will estimate the spin stiffness in the 
next section and corroborate this conclusion. Figure 5 
shows this behavior of Asp explicitly for the same pa- 
rameters as in Figure 4. 
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FIG. 5. Spin-reversed gap Asr in units of Ec as a function 
of T. The dashed line is Ez- Beyond about GK Asr is 
essentially the same as Ez- 

Now consider a noninteracting theory of CFs. This 
would have a completely flat dispersion for the collective 
mode, that is, w{q) = Ez- Therefore, as the temper- 
ature increases, the theory appears to be more weakly 
interacting, and the CFHF approximation becomes more 
accurate. This trend can be expected to continue until 
a temperature scale when CFs cease to exist. There are 
no obvious signs of such a scale in the data. 

The CF-Hamiltonian theory and the CFHF approxi- 
mation are very general, and can be applied to any frac- 
tional Hall state (for an application to compressible states 
see ref[42]). To ilustrate this Figure 6 shows the P{T) 
curves for v = 2/b for A = 1.5/ for a range of Zeeman cou- 
plings. Note the nonmonotonicity of the curves that start 
from the singlet ground state at T = 0. There is a tran- 
sition to the fully polarized state around Ez = Q-QlEc- 



been proposedEi.to explain intriguing observations by 
Kukushkin et ao of a state with half the maximal po- 
larization for v = 2/5, which is not allowed as a transla- 
tionally invariant CF state at T = 0. 

Let us compare our results to the only other method 
that can compute P{T) for arbitrary fractions, which is 
exact diagonalization (keeping all the excited states) aiid 
subsequent calculation of thermodynamic quantitiesCj. 
Due to computational limitations, this method is re- 
stricted to fairly small systems. For example, the latgest 
system studied by Chakraborty and PictilaincncS for 
V — 1/3 has 5 electrons, and for v — 2/5 has 4 elec- 
trons. For 1/ = 1/3 the exact diagonalization result lies 
above our predictions (and the experiment) for T > AK. 
This discrepancy might be the result of finite thickness 
and/or finite size corrections. However, at low T the ex- 
act diagonalization resulto follows the data more closely 
than our HF approximation (in all the above compar- 
isons the g — 0.5 line in Figure 2 of ref[48] has been used 
and compared to the lOW sample of Khandelwal et a£il. 
This sample has the closest parameters to those used in 
ref[48]). For v = 2/5, the CFHF results reproduce the 
nonmonotonicity of P(T) for those values of Ez where 
the singlet state is the ground state, and the peaks in 
P{T) occur at roughly the same^T in the CFHF and 
the exact diagonalization resultsa. However, the same 
overall pattern holds for i> =j-2/5, namely, the results of 
Chakraborty and Pietilainerc3 are below the CFHF pre- 
diction for low T, but are higher for T > 0.02Ec, where 
they once again see a 1/T tail with a large coefficient. It 
would be interesting to explore the finite size systematics 
to see if the large T tail is suppressed for larger sizes. 

IV. SPIN FLUCTUATIONS: MAPPING TO THE 
CONTINUUM QUANTUM FERROMAGNET 
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FIG. 6. Polarization versus T for v = 2/5 and A = 1.5?, for 
a range of values of Ez- Ez and T are both plotted in units 
of Ec- 

Note also that only results for translationally invari- 
ant HF states are presented, which ignores possible par- 
tially polarized inhomogeneous states that have recently 



As can be seen in Figures 1 and 2, spin fluctuations 
are quite important for v = 1/3 at intermediate tempera- 
tures, for ly = 1, the coarse-grained effective theory of the 
continuum quantum ferromagnet (CQFM) coupled with 
the large- approxirriation, first applied to this problem 
by Read and SachdevliJ, has been quite successful in ex- 
plaining the temperature dependence of the polarization. 
In this section a method is presented to map the = 1/3 
problem to the CQFM. 

In the CQFM description one starts with the actiontll 

l/T 

5 = Jd'^x J dT{iMoA{n) ■ V^n + ^{W^n)^ 



- MoH ■n+---) (14) 

where n is a local vector of unit length pointing in the 
direction of the magnetization, A{n) is the field that im- 
plements the Berry's phase needed to obtain the correct 
quantum commutation relations between the spin compo- 
nents, and H = g* ^bB is the Zeeman field {\H\ = Ez)- 
The two crucial parameters which enter the action are 
Mo, the magnetization density, and ps, the spin stiff- 
ness. There are other omitted terms in the action, which 
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produce at most a logarithmic correction to the mag- 
netization in two dimensions . There are various ways 
one can proceed at this pointer, but the most convenient 
one for our purposes involves mapping, the spin problem 
to a problem with Schwingcr bosonsE2l, and making the 
large-A'' approximationE2l. There are two common ways 
of mapping the spin problem to Schwingepiosons: The 
SU{N) approach and the 0{N) approachllJ, which give 
slightly different answers. 

In the following, we will restrict the theoretical results 
to the leading large-A^ approximation, where the magne- 
tization is given as M{T) — Afo$M(?', h), where r = Ps/T 
and h = Ez/T are scaling variables, and is a scaling 
functionj-Jn this approximation, the results of Read and 
SachdevEj for the magnetization in the SU{N) approach 
are 

^ f log(gi - e-^/^) - log(gi - e"/^) 

$Af (r, h) - — (15) 

and gi > 1 is a root of the equation 

(<Zi-e-'^/2)(<Zi-e'^/2) = g?e-«- (16) 



The corresponding results in the leading 0(N) 
approximatiortj are 



$M(r, h) = log(g^ " e-") - l"g(g2 ^ e^) ^^^^ 



where g2 > 1 is the solution of 

(g2-e-'')(g2-l)('?2-e'')=qfe-4- 



(18) 



The principal assumption underlying the CQFM de- 
scription is that long-wavelength ferromagnetic spin fluc- 
tuations are the only low energy modes that affect the 
spin polarization in the temperature range of interest. In 
the regime where this assumption holds for the v ~ \/2> 
state, it should be possible to map the CF-Hamiltonian 
theory to the CQFM. Conceptually, one can think of "in- 
tegrating out" the fermions and leaving behind an effec- 
tive theory of the spin fluctuations. Operationally, one 
needs to find the values of Afo and ps corresponding to the 
V = l/Z state. An additional complication arises here: 
Since the underlying fermionic theory responds to tem- 
perature by self-consistently modifying occupations and 
energies, one should expect to obtain temperature de- 
pendent values Mo(r) and Ps{T). We will extract these 
values from the CFHF Hamiltonian results. 

First consider Mq{T). We already have a value for 
the CFHF magnetization Mhf{T)- Let us first write 
down the correct relation between Mq and Mh f and then 
justify it. 

Mhf{T) = Mo(r)$M(r = 0,h = AsbJT) (19) 

In the CFHF theory the particles and holes are treated 
as independent, or noninteracting, with a gap equal to 
^SR (to be more precise, this is the lowest energy spin- 
flip excitation). To put it differently, the energy gap to 
create this spin-flip excitation is always Ag^^, no matter 



what the distance between the particle and the corre- 
sponding hole. This is as though the collective mode dis- 
persion were completely flat, Lj{q) — Asr- The CQFM 
description that corresponds most closely to the CFHF is 
the one that has the same spin-flip excitation spectrum, 
namely one with no spin stiffness, that is, r = 0, and and 
an effective Zeeman field E"/^ = Asr. The CQFM pre- 
diction for the magnetization of such a theory is the right 
hand side of equation (p^), which has to be equated to 
the CFHF prediction, hence the above formula, equation 

(ID- 
Knowing that, for SU{N), 



<i>M{0,h) = tanh(V2) 
while, for 0(iV), 

sinh(/i) 



$A/(0,/l) 



l/2 + cosh(/i) 



(20) 



(21) 



one extracts Mq{T) from equation ji^ , using the value 
of Mhp{T) computed in Section IIIEZT 

Next consider the spin stiffness ps (T) . At a given tem- 
perature T the self-consistent occupations NF,Gs{o'T'n) 
and energies e{a, n) in the ground state are computed 
using the procedure described in Section II. Now one cre- 
ates a twisted spin state by defining 

dhn.x = cos{qX/2)da.ji,x +sm{qX/2)d(3^n,x 

d],n,x = cos{qX/2)da,n,x + s\n{qX /2)di3^n,X (22) 

where a, (3 define local directions of up and down. In the 
twisted spin state the occupations of the local up and 
down spins remain the same as in the ground state, that 
is 



< d 



,X 



I >— 5nn'5xX'NF.GS 



(T,n) 



< d\,n,xdl3,n\X' >= ^nn' & X X' N F,Gs{i, n) 

< dl.n.xd/l.n'^X' >= 



(23) 



This leads to the following expectation values for the 
actual (global) t and | spin directions. 



< d 



d^ 



,X' > = 



Snn'SxX' X 



icos''iqX/2)NF,Gsit n) + sin^iqX /2)NF^GsiU n)) 



< d 



l,'n,X 



di,„',x' 



Snn'SxX' X 



{sin'iqX/2)NF,Gsitn) + C08^{qX/2)NF,GsiLn)) 

< d\,n.xdl,n' ,X' >= ^nn'^XX' X 

sm{qX/2) cos(gX/2)(A^F,Gs(T, n) - Nf.gs{L n)) 

< d\,n,xd],n',X' >= < d\ ,^ xdi.,n-,X' > (24) 

and corresponds to a state where the unit vector n point- 
ing in the direction of the local magnetization has compo- 
nents n = (sin(gX), 0, cos(gX)). Using the expectation 
values from equation (p4|), one can compute the HF en- 
ergy of the twisted ground state, and thence the excess 
energy to order q^. Comparing to the energy cost of a 
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twist in the CQFM, which is [ps/2)L'^q^ , one finds the 
spin stiffness 

{Npi^^n^) - 7VFa,ni))(iVf (T,n2) - NF{ln2)) (25) 

where is the area of the system, and Pmn^ is the den- 
sity matrix element of equation (|^). One caveat should 
be mentioned here: The above should be regarded as an 
estimate for the twist rather than a rigorous calculation 
(even in HF), since ideally one should compute the free 
energy cost of a twist, rather than just the internal en- 
ergy cost, as we have done. The free energy cost would be 
computed by carrying out a self-consistent HF at finite 
T in the presence of a twist. However, a fully relaxed 
HF calculation of an inhomogeneous state at T ^ is 
computationally prohibitive. In the following we use the 
value given by equation (|25|). 




FIG. 7. The parameters Mo{T) (normalized to 1 at T = 0) 
for SU(N) (dashed hne) and 0{N) (dot-dashed line), and 
ps{T) (x400 in units of Ec) extracted from the CFHF for 
Ez = 0.0186i5c and A — 1.51. While Mo shows no dramatic 
behavior, ps decreases precipitously beyond 3K. 

Figure 7 shows the results for Mq{T) and Ps{T) for 
the pacameters corresponding to the Melinte et al M242 
sampleE3. As can be seen, ps in particular has a dramatic 
T dependence, and essentially vanishes for T > QK. 
This corroborates the earlier observation that the spin- 
reversed gap collapses to Ez at around the same temper- 
ature, as seen in Figure 4. 

Having extracted the parameters Mo(T) and Ps{T), 
it is now easy to calculate the effects of spin fluctua- 
tions in the leading large- A'^ approximation from equa- 
tions ( p^ , p7[ ) . Figure 1 (in the Introduction) presents the 
results for the Khandelwal et al dataEJ. As can be seen, 
the agreement between the theory and experiment is now 
quite good, indicating that the proper treatment of long- 
wavelength thermal spin fluctuations has remedied most 
of the defects of the CFHF ncediction. Turning to the 
M242 sample of Mehnte et aO, a similar dramatic im- 
provement in the agreement between theory and exper- 
iment is seen in Figure 2 (in the Introduction). There 



are still some discrepancies in both the figures, but they 
never amount to more than 10 — 15%. 

One might wonder how important it is to keep the 
temperature dependence of the parameters Mq and ps. 
Figure 8 shows the prediction of the large- A^ approach 
using the values of AIq and ps computed at T = for 
the parameters corresponding to the Khandelwal et al 
dataliil. This fit ignores finite-temperature single parti- 
cle excitations. It is clear that while this prediction is 
somewhat better the the CFHF result for lower temper- 
atures, it is worse beyond 6K. Also, this prediction is 
uniformly worse than the one including both CFHF and 
spin fluctuation effects. 



1 r<m: 
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FIG.rSj The theoretical prediction for the Khandelwal et 
al datatij based on a zero temperature fit to the parame- 
ters of the CQFM. The thickness parameter is asswned to 
be A = 1.51. Circles represent the experimental datacil. 

Clearly, the reason is that this flt does not incorpo- 
rate the flnite-temperature fermionic single-particle ef- 
fects that renormalize the spin stiffness substantially 
downwards beyond about 3K, as seen in Figure 0. It 
must be concluded that the temperature dependence of 
the parameters Mq and p^ is crucial, and that both single- 
particle and thermal spin wave effects must be treated cor- 
rectly if one is to accurately predict the spin polarization. 



V. CONCLUSIONS, CAVEATS, AND OPEN 
QUESTIONS 

In summary, an approximate analytical method for 
computing the temperature dependence of the polar- 
ization for an arbitrary fractional quantum Hall state 
has been presented, and illustrated for = 1/3 and 
2/5. It consists of two steps. The first step is a finite- 
temperature Composite Fermion Hartree-Fock approxi- 
mation on the CF Hamiltonian. This already incorpo- 
rates important effects on the polarization resulting from 
particle-hole excitations, and is sufficient to produce all 

as the nonmono- 
For the spon- 



the correct qualitative features, sue 
tonicity of the polarization for 2/5 
taneously polarized 1/3 state, however, there is a sub- 
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stantial discrepancy between this prediction and experi- 
ment at intermediate temperatures due to the incorrect 
treatment of spin fluctuations. The second step, for the 
V = l/i state, is to map the long- wavelength low-eneijgy 
physics on to the continuum quantum ferromagneulj. 
The parameters that enter the effective theory are ex- 
tracted from the CFHF approximation. Finally, the ap- 
proximate solution to cpntinuum quantum ferromagnet 
in the large- TV approachlia is used to produce a prediction 
for the spin polarization which incorporates the effects 
of both single particle and spin wave excitations. Note 
that the reduction of the spin polarization due to single- 
particle and spin wave effects is not additive in the two 
effects: the parameters Mq and ps extracted from the 
CFHF approximation encode information about single- 
particle excitations in a very nonlinear way. The resulting 
prediction is in quite good agreement with experimental 
data, as seen in Figures 1 and 2. The comparison be- 
tween theory and experiment was made using the model 
potential of equation (|^) for illustrative purposes, with a 
thickness parameter A = 1.5Z, which ought to be in the 
physical regime for most samplesLJ. However, the above 
two-step procedure of first carrying out the CFHF, and 
then mapping the problem on to the CQFM, can be ex- 
ecuted for any potential v{q). 

A number of caveats must be noted at this point. In 
order for the mapping to the CQFM to make sense, spin 
wave excitations must be the only low energy spin-flip ex- 
citations. However, as T increases, other spin-flip excita- 
tions also become relevant (such as a transition from the 
f-spin n = CF-LL to the J,-spin n = \ CF-LL). Another 
complicating issue is that while the CQFM assumes the 
spin wave dispersion to be quadratic at all g, in reality it 
turns over and asymptotes at the spin-reversed gap Asfl. 
Both these considerations suggest that the mapping is 
to be trusted only in the regime where T < Asr{T). 
One can verify from Figure 5 that this corresponds to 
T < 6K or so, which fortunately includes the most inter- 
esting region of the data, ft is somewhat puzzling that 
the prediction from the mapping works well to twice this 
temperature. Finally, 1/N corrections can be expected 
to the leading large- iV prediction, which reduce the mag- 
netization slightly beyondjthe leading large- iV result (see, 
for example, Timm et atS). These corrections would im- 
prove the agreement between the theory and the data at 
low T. 

An important open problem is the development of a 
formalism in which one can systematically integrate out 
fermions, leaving behind a theory of the low-energy ex- 
citations, perhaps along the lines of Ref. [15]. Such ex- 
citations need not be restricted to spin wave excitations, 
but may perhaps include other spin-flip modes or even 
density modes, depending on the temperature and the 
system. Such a theory is necessary to address the effects 
of thermal spin fluctuations on the temperature depen- 
dent spin polarization of fractional Hall states which are 
not fully polarized in the absence of a Zeeman coupling. 
The spin-flip modes of such systems do not have the sim- 
ple quadratic CQFM form even at small q (for example, 
the spin- flip collective mode of the fully polarized 2/5 



state starts with a negative quadratic termed). A further 
interesting application would be to the interplay of spin 
and densibL. fluctuations in the compressible fractional 
Hall state^aa. 
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